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Mapping the Dirac point in gated bilayer graphene 
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We have performed low temperature scanning tunneling spectroscopy measurements on exfoliated bilayer 
graphene on Si02- By varying the back gate voltage we observed a linear shift of the Dirac point and an 
opening of a band gap due to the perpendicular electric field. In addition to observing a shift in the Dirac 
point, we also measured its spatial dependence using spatially resolved scanning tunneling spectroscopy. The 
spatial variation of the Dirac point was not correlated with topographic features and therefore we attribute 
its shift to random charged impurities. 



Monolayer graphene (MLG), which is just a single 
sheet of carbon atoms thick, has novel electronic proper- 
ties as a consequence of its linear band structure. Stack- 
ing one more layer on top of the monolayer gives rise 
to bilayer graphene (BLG) that is an exciting system 
with a different set of tunable properties^* 2 -. The bilayer 
structure is characterized by a quadratic dispersion re- 
lation E = ±h 2 k 2 /2m with the conduction band and 
valence bands touching making BLG a zero band gap 
semiconductor. When an electric field is applied per- 
pendicular to the plane of carbon atoms, it is possible 
to open up a band gap between the conduction band 
and valence band^£. Recent experiments with tech- 
niques like angle resolved photoemission spectroscopy 6 , 
infrared spectroscopy^^ and transport measurements 
with a double-gate^ have confirmed this band gap open- 
ing. These techniques are non-local and only provide in- 
formation about the average properties of the BLG. How- 
ever, from a device application perspective it is important 
to get details about how the spatial extent and morphol- 
ogy of the layers affect the electronic properties. Scan- 
ning tunneling microscopy (STM) is a powerful tool for 
this purpose. Previous STM studies have shown that im- 
purites in ML G 11 i 12 and phononsi&nfluence the charge- 
carrier scattering mechanisms in graphene. In this let- 
ter, we present scanning tunneling spectroscopy results 
for BLG on a SiC>2 substrate. These results show the 
spatial variation of the Dirac point as well as the control 
of the Dirac point and band gap due to the application 
of an electric field from the back gate. 

The BLG was prepared using the mechanical exfolia- 
tion techniqu e) 14 ! 15 . Degenerately doped Si with 300 nm 
thick Si02 on top was used as a back gate. Bilayer ar- 
eas were identified using an optical microscope and then 
Ti/Au electrodes were deposited using a shadow mask 
technique described elsewhere^. The device was then 
cooled to 4.6 K using an Omicron low temperature STM 
operating in ultrahigh vacuum (p < 10 -11 mbar). Elec- 
trochemically etched tungsten tips that exhibited a con- 
stant density of states on a Au surface were used for 



imaging and spectroscopy to avoid unwanted tip effects. 
Due to the cleaner fabrication procedure, no PMMA is 
used, it is possible to obtain atomic resolution images 
over large areas of the BLG without any additional clean- 
ing procedure unlike in previous STM measurements on 
exfoliated grapheneii*i2*i2*iL12*ia. 

Figure 1(a) illustrates the concept of band gap opening 
for BLG. The measurement set up with STM and back 
gate connections is shown in Fig. 1(b). The BLG is 
grounded and voltages are applied to the STM tip and 
the back gate. As the voltage is increased on the back 
gate or tip, an electrical field perpendicular to the BLG is 
established which shifts the Dirac point and opens a band 
gap. Figure 1(c) is a 60 nm x 60 nm STM image of the 
BLG showing modulations due to the underlying Si02 
substrate as typically seen in MLG samples on a SiC>2 



( a ) No electric field 
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FIG. 1. (color online) (a) Illustration of bilayer graphene 
and its band structure with and without a perpendicular elec- 
tric field. The electric field shifts the Dirac point and opens 
up a gap. (b) Schematic of mechanically exfoliated bilayer 
graphene flake showing gold electrodes on a SiCb substrate 
with STM and back gate connections, (c) Constant current 
STM image of bilayer graphene (0.5 V, 100 pA, 60 nm x 60 
nm). (d) Atomic resolution showing the triangular lattice (0.5 
V, 100 pA, 10 nm x 10 nm) 
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substrat o n i 12 i 13 i 17 i 18 i 19 . The triangular lattice arises due 
to the presence of two layers that leads to only a single 
sublattice being observed in the STM images 18 . 

Figure 2(a) shows a plot of dl/dV as a function of 
tip and gate voltage. The dl/dV point spectroscopy 
measurements are performed using standard lockin tech- 
niques with an ac modulation of 5 mV rms at 574 Hz. 
As we sweep through different gate voltages from -60 V 
to +60 V in steps of 2.5 V, two distinct observations can 
be made. The location of the minimum in the dl/dV 
curves shifts towards positive tip voltage with increasing 
gate voltage. Also, the width of the red region which 
corresponds to low values of dl/dV increases with an in- 
creasing gate voltage. In addition, there is an extra peak 
within the gap (white area visible at negative gate volt- 
ages) but further investigation is necessary to attribute it 
to any of theoretically predicted features like impurities 
and midgap states 20 . 

The gate voltage induces electrons or holes into BLG 
and generates an electric field perpendicular to the plane 
of the bilayer. It leads to an uneven induced charge on 
each layer and creates an electrostatic potential difference 
between the two layers that screens the external field and 
opens a band gap as predicted by the calculations^ 4 ^. 
The red region in Fig. 2(a) shows the band gap opening 
with gate voltage and the shift of the Dirac point. From 
the dl/dV curves we have determined the band gap at 
each gate voltage and the energy of the Dirac point based 
on the position of the minimum in the curve. 
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FIG. 2. (color online) (a) Point spectroscopy measurements 
of BLG as a function of back gate and tip voltage (energy). 
The color scale represents the local density of states, dl/dV, 
with red areas low and green areas high. The minimum in the 
density of states corresponds to the Dirac point and it shifts 
linearly with gate voltage. Similarly, the width of the min- 
imum increases with increasing gate voltage indicating that 
a band gap is opening, (b) Position of the Dirac point 
determined from the dl/dV curves in (a) as a function of the 
gate voltage V g (bottom) and induced carrier density n (top). 
A straight line fit to the data is shown, (c) Band gap varia- 
tion with V g . Dotted line is the fit (see text) to the band gap 
variation. 



A plot of the energy of the Dirac point as a function of 
the gate voltage is shown in Fig. 2(b). The variation of 
Dirac point with gate voltage is linear with a slope of 1.56 
meV/V. This linear dependence on gate voltage is differ- 
ent from MLG which has a square root dependence due 
to the linear bandstructure^. The linear dependence is 
expected for BLG because the bandstructure is quadratic 
E = h 2 k 2 /2m* where m* is the effective mass of electrons 
in BLG. For a two-dimensional system, the electron den- 
sity is written as n = k 2 /ir. Since the back gate and 
the BLG act as a parallel plate capacitor, the induced 



charge on the BLG is given by n 
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where t OT is 



the thickness of the oxide layer. Putting these equations 
together gives the shift in energy due to the back gate as 
E = h 2 7re€oV g /2m*et ox . A second effect of the gate volt- 
age is that it leads to an electric field between the two 
layers of graphene. Assuming no screening, i.e each layer 
carries charge density^ n/2, this field creates a potential 
difference between the two layers, V = | , 9 where d is 

ox 

the interlayer separation, 0.335 nm£ and the dielectric 
constant between the layers is 1 . This potential differ- 
ence reduces the effective voltage between the tip and the 
BLG. The net effect is that the energy of the Dirac point 
as observed by the tip shifts as 
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Using the above equation we find that the value of ef- 
fective mass that fits our observed shift of 1.56 meV/V 
is 0.023m e where m e is the mass of the electron. This 
value of the effective mass is slightly lower than tight- 
binding calculations which give m* = 0.033m e , however 
many body effects can reduce this mass^i. Furthermore, 
our analysis did not take into account the electric field 
due to the tip which would tend to lessen the effect of the 
gate and therefore increase the effective mass needed to 
fit our data. We have also neglected screening between 
the layers which reduces the potential difference between 
the layers^!. These effects could explain the discrepancy 
between our value and the calculated values of the effec- 
tive mass. 

Figure 2(c) shows the band gap as a function of gate 
voltage and carrier density n. The band gap has a 
finite minimum value around a gate voltage of 15 V. This 
is the voltage needed to have the Dirac point at the Fermi 
energy and therefore the minimum charge on the BLG. 
The magnitude of the band gap is fit using the equation 4 
A g = \e 2 VH\l(t\ + e 2 V 2 )\ x l 2 where V is the potential 
between the layers and the only free parameter is t± the 
interplane hopping. The best fit occurs for a value of 
t± = 0.12 eV which gives an effective mass of 0.017m e . 
Once again we have not taken into account the screening 
between layers which will reduce the potential difference 
and thus give a larger effective mass for our fit. 

In the case of MLG the shift of the Dirac point as a 
function of position leads to the formation of electron 
and hole puddles that can be seen in density of states 
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FIG. 3. (color online) (a) A 2D map (80 nm x 80 nm) show- 
ing the spatial variation of the Dirac point calculated from 
a dl/dV map of bilayer graphene. The colors correspond to 
positive and negative shift of the Dirac point, (b) dl/dV curve 
from the area of positive shift shown in (a) by a filled square, 
(c) dl/dV curve from the area of negative shift shown in (a) by 
a filled triangle. (d) Topography of the area shown in (a), (0.5 
V, 100 pA, 80 nm x 80 nm), the color scale shows the height 
variation of 1 nm. 

measurement :: 11 ! 12 ! 22 . To measure the spatial variation of 
the Dirac point for BLG, we have recorded LDOS maps 
for large areas of the flake, typically 80 nm x 80 nm. 
From the individual spectroscopy curves in the LDOS 
map we have calculated the energy of the Dirac point. 
This gives the energy of the Dirac point as a function of 
position in the BLG. Figure 3(a) shows the Dirac point 
for one such LDOS map. The color scale corresponds to 
the energy of the Dirac point. A dl/dV curve from an 
area of positive shift (red square in Fig. 3(a)) is shown 
in Fig. 3(b). Figure 3(c) shows a dl/dV curve from an 
area of negative shift (red triangle in Fig. 3(a)). The 
energy of the Dirac point is determined by the location 
of the minimum in these curves. The energy of the Dirac 
point is not evenly centered around Ej as the average 
value is at Ed = 20 meV. Figure 3(d) shows the topo- 
graph corresponding to the region in Fig. 3(a). The 
color scale here corresponds to the height variation of the 
image. We can clearly see from Fig. 3(a) and 3(d) that 
there is no correlation between the spatial variation of the 
Dirac point and the features in the topograph. Hence we 
claim that the origin of the Dirac point shift, at zero V g 
when no charge carriers are induced, lies in the random 
charged impurities as seen in the case of MLG. The impu- 
rity density variation can be calculated from the average 
measured shift in the Dirac point Ed using n = 2 ™ h 2 d 
to be 3.8 x 10 11 cm -2 . These results are consistent with 
an electrical transport measurement on the BLG which 
shows the Dirac point shifted slightly away from V g = 0. 
However, electrical transport measurements give only one 



number for the Dirac point and do not give any informa- 
tion about its spatial variation. 

In conclusion, we have presented scanning tunneling 
spectroscopy measurements on bilayer graphene on Si02 
at 4.6 K. Our measurements show that with a back gate 
it is possible to tune the band gap and doping in the 
bilayer. The opening of the band gap is a promising 
effect for tailoring the bilayer for electronic applications. 
In addition, the shift of the Dirac point due to random 
impurities emphasizes the need to advance the sample 
preparation methods to ensure minimum impurities. 
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